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Abstract 

I report the analysis performed on fermion masses and mixing, including neutrino mixing, within the 
context of a model with hierarchical radiative mass generation mechanism for light charged fermions, 
mediated by exotic scalar particles at one and two loops, respectively, meanwhile the neutrinos get Ma- 
jorana mass terms at tree level through the Yukawa couplings with two SU(2)l Higgs triplets. All 
the resulting mass matrices in the model, for the u, d, and e fermion charged sectors, the neutri- 
nos and the exotic scalar particles, are diagonalizcd in exact analytical form. Quantitative analysis 
shows that this model is successful to accommodate the hierarchical spectrum of masses and mixing 
in the quark sector as well as the charged lepton masses. The lepton mixing matrix, Vpmns, is writ- 
ten completely in terms of the neutrino masses mi,m2, and 7713. Large lepton mixing for #12 and #23 
is predicted in the range of values 0.7 < sin 2 20 12 < 0.7772 and 0.87 < sin 2 26> 23 < 0.9023 by using 
0.033 < s 2 3 < 0.04. These values for lepton mixing are consistent with 3a allowed ranges provided 
by recent global analysis of neutrino data oscillation. From Am 2 ol bounds, neutrino masses are pre- 
dicted in the range of values mi » (1.706 - 2.494) x lO^eF, m 2 » (6.675 - 12.56) x 10 -3 eV, and 
7TI3 « (1.215 — 2.188) x 10~ 2 eV, respectively. The above allowed lepton mixing leads to the quak-lcpton 
complementary relations 6g KM + 0[ 2 MNS » 41.543° - 44.066° and 0g KM + 9g 3 MNS « 36.835° - 38.295°. 
The new exotic scalar particles induce flavor changing neutral currents and contribute to lepton flavor 
violating processes such as E — > to radiative rare decays, r — > p/y,T — > e'y, (i — > ej, as well as 

to the anomalous magnetic moments of fermions. I give general analytical expressions for the branching 
ratios of these rare decays and for the anomalous magnetic moments for charged leptons. 
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1 Introduction 

The observed hierarchical spectrum of masses and mixing angles in the quark sector still remains as one of 
the most important challenges in particle physics. A possible solution to explain this hierarchical spectrum 
is that light fermion masses arise through radiative corrections pQ, while the masses for top quark, bottom 
quark, and tau lepton are generated either at the tree level like in Ref.[2], or by the implementation of 
seesaw-type mechanisms as was proposed by the author in a model with a SU(3) horizontal symmetry in 
Ref . [3J . Recently, it has also been possible to observe the phenomenon of flavor mixing in the leptonic sector 
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through the confirmation of the phenomenon of "neutrino oscillation" in experiments of neutrinos coming 
from atmospheric @j , solar [5] , reactors [6] and accelerators [7] , and as a consequence the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) lepton mixing matrix, Vpmns, has been determined. For three flavor neutrinos 
z/ e , z/,j, and u T , which may evolve into the three massive neutrinos v\ ,U2 ,^3, the experiments of neutrino 
oscillations are interpreted in terms of three mixing angles denoted by 612 for v e — v^ 823 for — i/ T , and #13 
for v e — v T . Recent analysis and fit of neutrino oscillation data give [8], at the 3a level, the allowed ranges 
of values 

I Am|, I = (1.4 - 3.3) x 10- 3 (ey) 2 , Am 2 2 = (7.1 - 8.9) x \^(eV) 2 , 

sin 2 26*23 = 0.87 - 1.0 , sin 2 29 l2 = 0.70 - 0.94 , (1) 

sin 2 6*13 < 0.051 , 

where Amf 2 and Am| are the solar and atmospheric mass differences, respectively. 

In this article I address the problem of fermion masses and mixing angles, including neutrino mixing, 
within the context of the model introduced in Ref.[2]. Section 2 briefly reviews the main features of the 
model with an U{1)h flavor symmetry. Next, in Sec. 3 I discuss the masses and mixing for charged leptons 
and quarks, at one and two loops, and using the strong hierarchy of masses, approximate mixing matrices 
for the u, d, and e charged sectors are provided. Section 4 is devoted to finding the upper bounds for mixing 
angles of charged leptons. In Sec. 5 I analyze neutrino mixing, and the Vpmns lepton mixing matrix is 
written in terms of neutrino masses, providing numerical results for neutrino mixing. In Sec. 6 I perform 
a quantitative analysis of quarks masses and mixing, including numerical values for the Vckm- In Sec. 
7 I give general expressions for the flavor changing neutral currents (FCNCs), rare decays and anomalous 
magnetic moments for charged leptons that are induced by the exotic scalar particles. Section 8 contents my 
conclusions. In the Appendix I have introduced a method to diagonalize in close analytical form a generic 
3x3 real and symmetric mass matrix, and then I have extended this method to diagonalize the 4x4 real 
symmetric exotic scalar mass matrix. 

2 Model with U(1)h flavor symmetry 

The gauge symmetry of the model is defined as U (1)h ® SU (3)c (8> SU(2) l ® f7(l)y. The fermionic content 
of the model is the same as in the "standard model" (SM), and their charges under the flavor symmetry 
U(1)h are arranged as to cancel anomalies without the introduction of exotic fermions. The fermions are 
classified, as in the SM, in five sectors / = q,u,d,l, and e, where q and I are the SU(2)l quark and lepton 
doublets, respectively, and u, d, and e are the singlets, in an obvious notation. 

The cancelation of anomalies in a simple way that simultaneously guarantees that only the third gene- 
ration of the charged fermions acquire masses at tree level is given by [2] 

H(f) = 0,±8 f , S 2 q -2S 2 u + 5 2 = 5 2 - 8 2 e , (2) 

with the constraints 

5t = 5 q = A^6 u = 6 d = 6 e = 6 (3) 

The assignment of flavor charges to the fermions is then as given in Table 1. The Gsm = SU(3)c <8> 
SU{2)l (g) U(1)y quantum numbers of fermions are the same as in the SM. 

The particle content of the model is such that we can implement a hierarchical mass generation mecha- 
nism, where the third family of charged fermions obtain mass at tree level, while the light charged fermions 
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Table 1: Assignment of family charges under U(1)h- 
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Table 2: Assignment of charges for scalar fields under U(l) H ® SU(3) C <8> SU(2) L <g> U(l) Y - 



get masses at one and two loops, respectively. In particular, the radiative mass generation for the light 
charged leptons involve the introduction of two SU(2)l weak scalar triplets with neutral fields, which we 
allow to get vacuum expectation values (VEVs). The VEVs of these triplets contribute very litle to the W 
and Z masses and simultaneously allow the generation of tree level Majorana mass terms for the left-handed 
neutrinos. 

The scalar fields introduced in the model are then divided into two classes. Class I ( II ) contains scalar 
fields which acquire (do not acquire) VEV. These scalar fields are as given in Table 2. 

The Yukawa couplings are classified in two types, Dirac(D) and Majorana(M) [Figs, (la) and (lb)], 
respectively, 

Cy = C-yd + -Cym , (4) 

where 

Cyd = Y u q L3 4> lUR3 + Y d q L3 fad R3 + Y e /" L 30iTR3 + H.c , (5) 
with * = ^7°, = io20*, Y\ where i = u, d, e are coupling constants, and 

Cym = Y\ 2 qfJCfo^q^ + Yhs C ' fa {a ,f3 } q§ L + Y« 33 q% L T Cfo^q^ 

+ Y u 12 u 1R T C(f) 7 U2R + Y u 23 u 2R T C(j)8U 3 R + Y u 33 u 3R T C(f) 7 u 3R (6) 
+ Y d 12 d 1R T Cfod 2R + Y d 23 d 2R T C(j )6 d 3R + Y d 33 d 3R T C^d 3R + h.c. 

+ Y u if L T C4> 10{at(3} i^ L + Y 23 i 2L T C4> 9{at(3} i^ L + Y 33 i 3L T C4> 10{at(3} i^ L ^ 

+ Y 12 e R T C(\)\ 2 [i R + Y 2 3 HR T C<t>uT R + Y 33 T R T C(f> 12 T R + h.c. 

In these couplings C represents the charge conjugation matrix, a and are weak isospin indices, and color 
indices have been omitted. Couplings of Eq.(6) are introduced for the quark sector, while those of Eq.(7) 
are needed for the lepton sector. 

Notice that 03 and 0g are represented as 



3 



-4/3 0-1/3 
-1/3 ^2/3 



and 



k++ 



(8) 



where the superscripts denote the electric charge of the fields ( and corresponding expressions for 04 and 
010 )• 

The most general scalar potential is written as 

-v(4>i) = ^2i4\\<t>i\\ 2 + ^2^ij\\<f>i\\ 2 \\4>j\\ 2 + V3i<t>l<f>i<f>3<t>i + mi4>\4>l<p3<t>i 

i i,j 
A$ A. A. _L A$ a3 A. 1, _L sh„A, _L Z-„, A$ A$ . 



+?7410i040401 + ?741 01040401 + K 9 i010^090i + K9101090901 

+Kio,i0i0io0io0i + Kio,i0i0io0io0i + ^ ^ij 110^0:? I| 2 

+ ( P1050602 + P2070802 + Ai0^0" 3 {a,/3} 01 + A20 7 0?0 3 {a,/3}01 
+A 3 Tr(0j0 4 )0| + A 4 050 6 7 02 + A 5 5 0j0 7 8 + A 6 2 080l + 2/1 0120H02 



+ Cl01201 01O{a,/3}01 + ^r7>(0{ o 9 )02 + ei0 5 0^0| 2 
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(9) 



+ €20708012011 + h.C. ) 

i^"0i an appropriate contraction of the SU(2)l and SU(3)c indices 



where Tr means trace and in ||0j|| 
is understood. The gauge invariance of this potential requires the relation A = 25 to be hold. 

In order to break the symmetry, the VEVs for the class I scalar fields are assumed to be in the form 
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, (02 ) = V 2 , 
, (010) = 



vio 




(10) 
(11) 



(01 ) and (02) achieve the symmetry breaking sequence 



U(1) H ® G SM ^ G SM ^ SU(3) C U(1) Q , 



(12) 



while the VEVs vg and v\q are extremely small in order to be consistent with the experimental bounds on 
the p parameter. My/ = \gv\ with v\ 246 GeV, and I assume v 2 in the TeV region. 

The scalar field mixing arises after spontaneous symmetry breaking (SSB) from the terms in the potential 
that couple two different class II fields to one class I field. After SSB the mass matrix for the scalar fields 
of charge +2 , (0g, 0io, 0i2, 0n) may be written as 
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(13) 



V>2 e{ x I 



2 

where e| = /xf + Aji-y + -^2^1 l° r ^ = 9, 10, 11, 12, and analogous ones for the — I and | scalar sectors. 
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3 Masses and mixing for charged fermions 



Now I give a brief description of the hierarchical mass generation mechanism for the charged fermions. After 
the SSB of the electroweak symmetry down to U(1)q of QED, the Yukawa couplings of Eq.(5) generate 
tree level masses for the top and bottom quarks and the r lepton. For the light charged fermions, the scalar 
fields introduced in the model allow the one and two loop diagrams of Fig. 2 for the charged lepton mass 
matrix elements, and similar ones for the up and down quark sectors. In the diagrams of Fig. 2 the cross in 
the internal fermion lines means tree level mixing and the black dot means one loop mixing. The diagrams 
of Figs. 3(a) and 3(b) should be added to the matrix elements (1,3) and (3,1), respectively. 

In the one loop contribution to the mass matrices for the different charged fermion sectors only the third 
family of fermions appears in the internal lines. This generate a rank 2 matrix, which once diagonalized gives 
the mass eigenstates at this approximation. Then using these mass eigenstates the next order contribution 
is computed, obtaining a matrix of rank 3. After the diagonalization of this last matrix the mass eigenvalues 
and eigenstates are obtained. 



3.1 Charged leptons 

The nonvanishing contributions from the one loop diagrams of Fig. 2 to the mass terms eiRejL^i/ 1 ^ are 

S 22 = m rJ^2 E U lk U 4k f(M k ,m° T ), (14) 
k 

^23 = m °^f T, U ^kf(M k ,m° T ), (15) 

k 

4? = ™°ll|r E UikU 3k f(M k , m° T ) , (16) 

k 

where m° is the tree level contribution, U is the orthogonal matrix which diagonalizes M|, with 

<f>i = UijVj , $1 = 09 , $2 = 010 , $3 = 012 , $4 = 011 (17) 

being the relation between gauge and mass scalar eigenfields cij, i,j = 1,2,3,4, M| are the scalar mass 
eigenvalues, and 

/M)-^ 1 ^- (18) 

Thus, the one loop contribution to fermion masses may be written as 

\ 

(19) 
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Now, Me is diagonalized by a biunitary transformation 



V^MWV™ = Mj» , (20) 
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where in this report (ignoring CP violation), the orthogonal matrices and V R L> are 





a R c 'r 
c' R b' R 



where 



(i) 



1 

cosa^ — sinak 
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(1) 



R 



1 

cos aR — sin or 
sin or cos aR 



cos«l = c/(A+ — a' L ) 
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X+-X- 
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A_ and A+ are the solutions of the equation 
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^23 ~ "32 



D> = a L b L ~ C> L = a R b R ~ C R 
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0203 - 023032 > , 
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VP T M^ T M^VP = V^ T mPm^ T V { r X) 





1 I assume the signs: 0,2 > 0, (223 < and 032 < 



Therefore, from Eq. (30), up to one loop level 

m W = , =yf\Z , m 3 1} =7a7 , (31) 

with the expected hierarchy A- <C A+. 
Two loop contributions for charged leptons: 

S n ^E^'^W^'l^WtM^f'), (32) 



4? = ^^^^rnf\v^\ i {V^\ i U lk U zk f{M k ,mf ) ), (33) 

k,i 



+^E M ^ 1) (^ 1) )3 i (^ 1) )2 i C/ 2fe C/3 fe /(M,,m; : 



(34) 



4? = E^ 1) ( y i 1) )^( V K 1) )3^2 fc f/ 4fc /(M fc ,mJ 1) ), (35) 

k,i 



4? = ^E M ^ (1) (^ 1) )2,(^ 1) )2^ 2 ^ 4fc /(M fc ,m; i) ) 

+^ Y. k ,r4\vi^(V^U2kU, k f(M kl m?) , 



(36) 



16vr 

where i = 2, 3 and fe = 1, 2, 3, 4. 

Hence the mass matrix for charged leptons up to two loop contributions may be written in good approx- 
imation as 



Mi 2 ) 
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a/a7 I = I S 21 y/A. I . (37) 
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2 / -I \ 

In the limit M k » the function f(a,b) behaves as hi^. In this limit, and introducing the m\ one 
loop mass eigenvalues, Eq.(31), the orthogonal matrices and V^, Eq.(23), the relationships 



V V nr\c 1 1/ r cm r\i ^ — — 1 1 

(38) 



cos oll cos = — a + -a — ' cos oil sin = — a + -a 



a2\/A+— as W\- . a32\/A++023\/A- 

sm «i sin otR = — v A+ _ A v — , sinaicosafi = v A+ _ A — - — 

and using the orthogonality of U, one obtains explicitly 
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E n = a 2 o > 



] l2 = E 2 1 = - ^3^32 > C = ^3 

1Z Zi c a a 3 ' Q Y12 



Si3 = c a a 23( 7+1^<0 , 
£3i = c a a 32 a + ^^f <0 , 



where the parameter a is defined as 

Y 2 M 2 

A; 

3.2 Quarks 

Performing an analogous analysis, the one loop contributions for the down quark sector are 

, d nYJLY£, 



(S«) =rn°^J2u? k Uif(M*,rn° b ) = a d 2 , 

k 



k 



fc 



where m° is the tree level contribution, U d is the orthogonal matrix which diagonalizes Aff ^, with 



being the relation between gauge and mass scalar eigenfields erf, i,j = 1,2,3,4 and Mf} 2 are the 
eigenvalues. 

Similarly, the two loop contributions for the down quark sector may be expressed as 



= a\ a d > , 
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r 33 m b 



(45) 
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*33 ^C? I ^12 a 2 a 32 ^ A 

1 12 *33 m fe 



with <r d defined as 



yd 
^31 



^33 _i_ ^12 Q 2 a 2l 



Mi 



(46) 



and analogous ones for the up quark sector. 

Notice that the radiative corrections give rise to the following hierarchy among the parameters in the 
mass matrices for charged leptons 



En , E12 , E 2 i , |Ei 3 | , |E 3 i| <C a 2 , |a 23 | , |a 32 | < a 3 = m° , (47) 

and similarly for the u and d quark sectors. 

The matrix M^, Eq. (37), is diagonalized by a new biunitary transformation 

V^M^V^ = M$ . (48) 

(2) 

In this case, with the aid of results of the Appendix, the orthogonal matrix V£ may be written as 



/ l{bL-\l){CL-\l)-fl 



V, 



(2) _ 



(A 2 -A 1 )(A 3 -A 1 ) 



(«L-Ai)(cL-Ai)-e| 
(A 2 -Ai)(A 3 -Ai) 
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\ V (A 2 -Ai)(A 3 -Ai) 
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(b L -x 3 )(cL-x 3 yfi \ 

(Ai-A 3 )(A 2 -A 3 ) 



>L-A 3 )(cL-A 3 )-e| > 
(Ai-A 3 )(A 2 -A 3 ) 



(aL-A 3 )(bL-A 3 )-4 
(A1-A3XA2-A3) / 



(49) 



where A« ,i = 1,2,3 are the eigenvalues of ( AfjP ); Ai = ml, A 2 = m^, and A3 = m 2 T 

(2) 

for charged leptons. A similar expression is obtained for by replacing L by R, where the L and R 
parameters are defined as 
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. a L d L e L \ I a R d R e R 

MP MP = [ d L b L f L , M^mP ^ d R b R f R ) . (50) 
eL Jl c L / \ e R f R c R 



a L = E?i + E 21 + So, , c?l = E 2 i\/AT + S11S12 , ej, — E 31 -v/X^ + E n E 



11 "r ^21 "r ^31 i «L — ^21 V A - "T ^11^12 5 e £ — ^31 V A + ^11^13 

6 L = A_ + Ef 2 , / L = E 12 S 13 , (51) 

c L = A+ + E? 3 

a R = £-11 + E^ 2 + E^ 3 , d R = E^a/aT + E11E21 , e R = E13 a/A^ + E11E31 

6 i? = A_ + E| 1 , / R = E 21 E 3 i , (52) 

c R = A+ + E§! 

and similarly for the u and d quark sectors. 

Up to now I have realized a complete exact analytical diagonalization of the resulting charged fermion 
mass matrices at one and two loops contributions. Thus, from Eqs. (23) and (49) one may obtain exact 
expressions for the orthogonal matrices 

V L = vjpvjn , Vr ^V^V^. (53) 

In particular, the Vckm quark mixing matrix (ignoring CP violation) takes the form 

Vckm = (V^fvM (54) 

3.3 Approximate mixing matrices for charged fermions 

Taking advantage of the strong hierarchy of masses observed in the quarks and charged leptons, it is possible 
to make good approximations for the mixing matrices and given in the Eqs. (23) and (49), 
respectively. For instance, the tree level contribution defines the magnitude of the masses for the heaviest 
fermions mt , mj, and m T in each sector. One loop contribution determines the masses for m c , m s , and m M 
and the mixing angle between the 2 and 3 families, giving simultaneously small corrections to the masses 
of the heaviest fermions. Finally, the two loop contribution gives masses to the lightest fermions u, d and 
e and determines their mixing with the families 2 and 3, giving some tiny corrections to the masses of the 
heavier fermions. I use this perturbative mechanism to make the following approximations: 



, . .. D' X-X, 
A_ <C A + implies 



A_ 



(A_+A + )2 A + (i + A_)2 
and expanding now the square root in Eq.(26), one gets 

~ n. 

~ B' 



«1, 



(55) 



A_ = f 



D' 



(a2a3-a23Q32) 2 



{a 2 



s 12 ) 2 



(56) 
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and hence from Eqs. (21) and (27), \ + — b' L = a! L — A_ « a| + a§ 2 — (02 — S12) 2 ~ a§ 2 . Thus, from Eqs. 
(23)-(27), the strong hierarchy of masses leads to the approximation 



V, 1 ^ ( D C 23 S 23 ) , s,, , . s /i -4,. (.-,7) 

S23 c 23 1 V A;i ^ 

where I have set A+ — A_ ~ A3 — A 2 . Then, in this approach the main contribution to the mixing angle between 
the families 2 and 3 comes from the mixing matrix element {V^)2Z- So, according to this perturbative 

approach the two loop contribution (V^ 2 "*)^ is negligible, that is, KV^^I <C KV^^I, and hence it is 

(2) 

a good approach to set (V£ J23 ~ in Eq.(49). Let me discuss further about the consistence of this 
approximation. 

From the results of the Appendix, by replacing a,b,c,d,e, f — > ai, &l> c l, di, e^, /l, one gets the exact 
equation 

A 2 l(A 3 ) = A3 - (a L + c L )A 3 + a L c L - e 2 L = {a L c L - e| - A1A3) - A 3 (A 2 - b L ) , (58) 



where a L c L -e\ = (Ef x + X 21 )A + - 2Eii£i3E 3 i + (S^ + S|i)S? 3 ~ 0(2 loops) 2 A+. Let me recall here 
that from Eq.(31), A 2 — > A_ and A3 — ► A+ if one ignores two loop contributions to the masses for the 2 and 
3 families. Therefore, because Ai ~ 0(2 loops) 2 and A 2 — = A 2 — A_ — £f 2 ~ 0(2 loops) 2 by construction, 
one concludes that the right-hand side of Eq. (58) is ~ 0(2 loops) 2 A3, and hence 



One can also directly write 



A| - (a L + c L )A 3 + a L c L - e 2 = A 2 - [A_ + A+ + a L + £ 2 3 ]A 3 + [A_ + ^-^]A 3 

« A|-[A_ + A+]A 3 + A_A 3 



(60) 



Comparing now Eq. (60) with Eq. (25), one concludes, in what concern the magnitudes for the heaviest 
fermion in each charged sector; m 2 , m 2 , and m 2 , that the eigenvalue A3 satisfies with good approximation 
the quadratic equation 

A 2L (A 3 ) = A| - (a L + c L )A 3 + a L c L - e\ » , (61) 

and hence (V£ ) 2 3 ~ is in good agreement with the radiative corrections. Using now Eq. (61) and the 
orthogonality of vj® , Eq. (49) , one reaches the approximation 



, Ci 2 Ci 3 Ci 3 Si 2 -S13 \ / c 13 -S13 \ / C12 S12 

^ 2) « I - Sl2 Cl2 0=010 -s 12 cia I , (62) 

Cl 2 Si 3 Si 2 Si 3 C13 / V s 13 C13 / \ 1 



where the si 2 and S13 mixing angles are identified as 
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^21 / ^31 

S 12 = \l X^xl > S 13 = A/ aPa^ 



(63) 



Cl2 = a/1 - S? 2 , C13 = a/1 - S? 3 , 

and therefore 



C12C13 C13S12 — S13 

L2 — C12S13S23 C12C23 — S12S13S23 — C13S2; 

C12C23S13 — S12S23 C12S23 + C23S12S13 C13C23 



Vl = Vtpvj® ~ I — C23S 12 — C12S13S23 C12C23 — S12S13S23 — C13S23 I • (64) 



It is important to emphasize here that Eq. (64) it is not a parametrization but an approximation that 
is consistent with 

• the strong hierarchy of masses for quarks and charged leptons and 

• the radiative corrections. 

A similar analysis for the R handed mixing matrices yields 

■>' ' I I I II II / / / / 3 

VR = Vr Vr ~ — C 23 S 12 — C 12 S 13 S 23 c 12 c 23 — S 12 S 13' S 23 — C 13 S 23 i (65) 



I s~J fj s~J I s~J /~'^ 

5 12*23 c 12 c 23 5 13 c 12 5 23 ' L 23 5 12' b 13 C 13 L 23 / R 



where 



■ 'WA ■ s - s vA- (66) 

So, the strong hierarchy of masses for charged leptons leads to the approach 
Sl9 ~ Si ~ 2loo P s ^l s ' « ^ « ( 2loo P s ^l 

12 ~ m M ~ ^ x l 00 p y ' 12 ~ m M ~ ^ x l 00 p / ' 

S9 „ ~ W w ( iloop \ / w k«l w ( iloop \ ( 67 ) 
23 ~ m T ~ \tree level / ' 23 ~ m r ~ \ tree level / ' 

s ~ ISail w Q ( 2 loops \ , ^ [Siaj ^ q / 2 loops \ 

13 ~ m T ~ V tree level/ ' 13 ~ "v ~ \ tree level/ 
Hence one concludes that the radiative mass mechanism naturally yields the hierarchy 

S13 < S12 , S23 S S'l 3 < s 'l2 , s 23 • ( 68 ) 

Notice that the approximations realized in this section, Eqs. (57) and (62)-(68), may equally be applied to 

the e,u, and d charged sectors by replacing properly the mass parameters: — ► , m° , 02, 023, 032 — ► 
a 2' i a 23 j a 32 an< i ^11) ^12) ^13) ^21, ^31 , , ^13 , ^21 > ^31 > respectively, in an obvious notation. 
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and hence from Eqs. (69)-(72): 



a 2 



4 Quantitative analysis of masses and mixing for charged leptons 

Using the orthogonality conditions of U and Eq.(171) of the Appendix, one can compute the parameters 
involved in the charged lepton mass matrices; one gets explicitly 

«2 - m rl 67 r2 Z^fc=l u lk^4k Ul — "V 167r 2 Z-,fc=l 111 — "V 16tt 2 ^ 22 

, (69) 

„ _ ^0 ^23^33 V 4 TJ,TJ,\-n M k — m° Y23X33 V 4 /sfafc) l n % — m 0Mk /? 

_ y i 2 2 V 4 77,77, In M * - V 4 33(^) ln % = yi22 F 

where % = M|, and the dimensionless functions F22, -F23, -F32, and Fo- may be written as 
^22 = Yi ^Mln^ = ln ^ - ln 2* + ln 2± > 

F 23 = - V? „ 9 -^ln^ = ^Fxr cd t~? ] + ^r^r In 21 > 

F 32 = _£ 4 M^/ n 2* = tr?i' ^ ^ - b A^zml ln 33 + 6c(°4-»?4) ln 2i > ( 7 o) 

p _ V^ 4 93iVk)_l n m _ c(a±-ri2){a4.-m) J n V2 _ c(ai-r?3)(a4-q3) j R »?3 

c Z^fc=2 h(»? fe ) ?7i IM»?2)| r/i /i(r/ 3 ) 771 

I c(ai-q 4 )(a4-q4) j n q 



where 

\Km)\ = im - m){m - m){m - m) , \Km)\ = {m-m){m-m){m-m)- (71) 

To leading order in the radiative loop corrections, one reaches the approximations 

m T = V^3 ~ a/A+ ~ % = 1 

"V = V 7 ^ ~ ~ «2 , ( 72 ) 



-^22 ~ O — TFZ?, fir ~ ~ 

or (73) 

^ -243842 . Fli ^ .0695437 
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So, in this approach the following relationships hold: 



112*23 V ^7 _ 2.875643 x 1Q- 4 
16tt 2 J ~ F 22 F a F 22 F a 



(74) 



^ « (.285199)^ , (75) 



|<3.23 | 


~ |<*23| 


^33 ^23 


r ^12-^23 






^23-^22 


Ca Y 23 F 22 


|«32| 


~ l«32| 


^33^32 


£ *12 -F32 




02 


^23-^22 


Q 5^32 F 22 



J23 



^3: 



(76) 



v, 1 «23fl32 1 |«23 1 1 032 1 1 1 1 e 1 , , /e 

2^12 — 2^21 — — q ~ — 1 0-23 1 ^23 = — l a 32| «23 1 {'') 

C a 0-3 C a TTl T C a C a 



1 2.31 1 - c a \a 3 2\<J + -— c a — |a 32 | + -^l a 2 3 | 



|S 13 | = c a |a 23 k + 12^ « c a ^|a 23 | + ^|a 32 | 



(78) 



where the superscript e denotes the charged lepton sector. Therefore, the mixing angles in Eq. (64), 

and (V e R), Eq. (65), may be expressed as 



4 3 « \/5 «^ = ^ *b = c a (0.016957)^ 

S ? 2 - ^/^«^ = ^*« (-285199)^^ 



(79) 



S 13 ~ VA,-Ai ~ m. ~ C "m„ S 23+UJ ^ 



1 2:s 



and 



s i3 ~ v x^xr ~ ~ Ca m7 s 23 + v^TTT/ z: 



A 3 — Ai m r a m„ ^ 



'23 
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respectively, with the relations 



'12 



= S 



le 
12 



and hence 



23 
/ 

23 



|«32| 
I «23 1 



^32 
^23 



(81) 



4.1 Upper bounds for charged lepton mixing angles 



Each particular set of scalar mass parameters in Eq.(165) define a spectrum of scalar mass eigenvalues r/i, 
V2 , f]3i Vii the values for F22, F23, F32 and F a through the Eq.(70), as well as the magnitudes for mixing 
angles in V e L and V e R through Eqs. (79) and (80). A numerical evaluation shows that the variation of these 
mixing angles is relatively small for a large region in the space mass parameters. So, in order to find out the 
orders of magnitude for these mixing angles, let me redefine the parameters of M?, Eq.(165), in such a way 



Mi 



b' 



\ 


d! 



(82) 



\ d' a' A J 



where a[, a' 2 , a' 3 , a' 4 , b' , d and d' are positive real numbers, while M is a mass parameter in the TeV region. 
Mixing angles in V e L and V e R do not depend on M. The value of M 2 may be determined for instance by 
specifying the value of the lightest scalar mass eigenvalue 771 = Mj\ 



Setting for example and simplicity a[ 



V 



d' one gets F23 = -F32 and hence s 



23 



Je 
s 23> 



6 12 



s l2i s 



13 = s 'i3- i n the simplified parameter space region defined by a[ = a' 4 = 1, 2 ^ a' 2 S 120, 
3 ^ a' 3 ^ 125, 1 ^ V = d = d! ^ 10, one gets the following range of values for mixing angles in the 
charged lepton sector: 



4.666163 x 1(T 3 < 


5 23 


Je 

— *23 


< 1.34417 x 10 _ 


0.105688 < 


„e 
*12 
„e 
A 23 


„/e 

s 12 

„/e 
*23 


< 0.240805 , 


4.931603 x 10~ 4 < 


S 12 


— s /e 


< 3.236832 x 10 


3.904076 x 10~ 4 < 


S 13 


Je 

— 6 13 


< 9.123709 x 10 



(83) 



where the upper and lower bounds are obtained with the values a' 2 = 2, a' 3 = 3, b' = d = d! = 1, 
and a' 2 = 120, a' z = 125, b' = d = d! = 10, respectively, and where I have used the range of values 
° a = T?f = 2I1 = 0-742528 — 0.938792 corresponding to the global parameter space region defined by Eq. 
(103) in the analysis of neutrino mixing. 
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5 Neutrino masses and Vpmns 



From the Yukawa couplings of Eq.(7), the mass matrix for the left-handed neutrinos is obtained as 

/ Y 12 vio/2 \ / d \ 
M v = Y 12 v w /2 Y 23 v 9 /2 = id / . (84) 
\ Y 23 v 9 /2 Y 33 v 10 J \ / c / 

One may diagonalize M v as 

UjM u U u = Mt , (85) 

where = diag(^i,^ 2 ,^ 3 ) is the diagonal matrix with £i , £ 2 , and £3 being the eigenvalues of My, and C/^ 
is the rotation matrix which connects the gauge states with the corresponding eigenstates. 

The eigenvalues £1, £ 2 , and £3 satisfy the following nonlinear relationships with the parameters d, f and 
c of M v , Eq.(84): 

6 + 6 + £s = c 
£16 + 66 + 66 = ~d 2 -f 2 (86) 
666 = -d 2 c 

The square matrix elements Uyfj may be obtained from those of (V^ )fj-, Eq.(49), by replacing 
a<L,bL,CL,dL,eL, Jl — > 0,0, c, d, 0, / and Aj — > & respectively. However, from the Eq.(86) and assuming 
c,d,f > 0, it is easy to conclude that one of the 6 z = 1,2,3 is negative. Thus, the eigenvalues 6 cannot 
be directly associated to the physical neutrino masses. 

Setting £3 > and computing explicitly the Ujjj elements, one arrives to the following statements: 

1. Assuming normal hierarchy 

£1 < £2 < £3 implies £1 > and £2 < , (87) 

2. Assuming the hierarchy 

£2 < £1 < £3 implies £1 < and 6 > 



In what follows, I assume a normal hierarchy for the squared eigenvalues £f as in Eq.(87@. In the 
literature there exists a lot of models dealing with normal neutrino mass hierarchy [9]. Now I define the 
unitary matrix V u = U u diag(l,i, l)|10j. which may be written as 



Vy 



P m 2 m 3 



(m|-m?)(m§- 



f 2 c mi 
(m\— mf)(m^— m'f) 



{(P-mpjdP+p-ml) 



P mim z 



(m^—mylm^—m 2 ,) 



P c 777,2 



{m 2 -d 2 ){d 2 +p-m 2 ) 
(m 2 -m 2 )(m 2 -m 2 ) 



P 7711777-2 



(m 2 —m 2 )(m 2 —m 2 ) 



f cm 3 
(m^— mf) (777,3 — m 2 ,) 



(m 2 ~d 2 )(m 2 -d 2 -p) 
(?7j|— m 2 )(m 2 — m|) / 



The second possibility, Eq. 



does not change any conclusion about this model 
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or equivalently in the form 



m 2 m 3 {m 3 -m 2 ) 
c(m 1 +m 2 )(m 3 — mi) 



TOi(m 3 — m 2 ) 
(mi+m 2 )(m3-mi) 



TO 1 (m 1 +r7i 3 )(m2— mi) 
c(m 1 +m 2 )(m 3 -TOi) 



miTO 3 (mi+m 3 ) 
c(TOi+m 2 )(m 3 +m 2 ) 



m 2 (rrii+m 3 ) 
'm 1 +m 2 )(m 3 +m 2 ) 



m\m 2 {m 2 —m\) \ 
c(TO 3 -mi)(m 3 +m 2 ) 



m 3 (m 2 -mi) 
(m 3 -TO 1 )(m 3 +m 2 ) 



m 2 {m 3 -m 2 )('m2-m 1 ) / m 3 {m 1 +m 3 ){m 3 -m 2 ) 



c(mi+m 2 )(m3+m 2 ) y c(m 3 -rni)(m3+m 2 ) / 



(90) 



after using properly the results of the Appendix and making the identification 

(6b - 6,6) = (roi,m 2 ,m 3 ) 



(91) 



between the eigenvalues £j and the physical neutrino masses mi, m2, and 7713. Therefore, for neutrinos the 
transformation between gauge, = (v®,v®,v®)l, and mass eigenstates, V^l = (^b ^2, ^3)l 5 is 

V>° L = K ^ . (92) 

From Eq.(91) and the definition of V v it is easy to verify that 

V„M V V V = diag(mi,m 2 ,m 3 ) , (93) 
and hence one may write Eq.(86) in terms of neutrino masses: 



mi + 771,3 — m 2 
m\m,2 — mirri3 + m^raj, 
mim2Tri3 



d 2 c. 



(94) 



The combination of these relationships yields the useful equality 

f 2 c = (m 3 - m 2 )(m 3 + mi)(m 2 - mi) 



(95) 



Notice that Eqs. (94) and (95) allow one to write all the matrix elements (V u )ij completely in terms of 
the physical neutrino masses mi ,771,2 , and 7113 as in Eq.(90). 



5.1 Vpmns lepton mixing matrix 

The current experimental study of neutrino oscillation phenomena gives as a result that in the lepton sector 
the mixing matrix Vpmns behaves close to the so-called "tribimaximal mixing" (TBM)[llJ. In particular, 
according to Eq.(l), the mixing angles 9\2 and 623 are large, sin #12 and sin #23 < O(l), while #13 has not 
yet been measured. So, taking the ranges of values in Eq.(83) as the typical orders of magnitude for the 
mixing angles in the charged lepton sector, it is then clear that mixing in the lepton sector should come 
almost completely from neutrino mixing, and then one may approach with good precision 



Vpmns = {V eL )W v « V v . (96) 
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Thus, from Eqs. (89), (90), and (96), the Vpmns lepton mixing matrix in this model may be approached as 



C12C13 -% C13S12 S13 

Vpmns ~ I C23S12 - C12S23S13 i (C12C23 + S12S23S13) C13S23 
- (S12S23 + C12C23S13) ~i (C12S23 - C23S12S13) C13C23 



(97) 



where the lepton mixing angles are identified as 



^13 — (K)i3 — 



f 2 - 

(m 2 -m 2 )(m|-m2) 



f 2 TO1TO2 



mim 2 (m 2 — m-i) 



(TOi+m 3 -TO 2 )(m 3 -TOi)(TO 3 +m 2 ) 



g2 — m\-m\ p mi m 3 

12 m 2 (m 2 -m 2 )(m§-m^)-/ 2 mim 2 



_ sf 3 m 3 m\-m\ 
c 2 3 m 2 m\—m\ ' 



(98) 



^23 — 



/ 2 cm 3 



(m|— m 2 )(TO 2 — to 2 )— / 2 TO1TO2 



£13 EM _c_ 

c 2 3 TO 2 TOl ' 



and 



^3=1 



2 (to 2 — to 2 )(to 3 — to|)— f 2 mim 2 

; 13 - (to 2 -to 2 )(to 2 -to 2 ) ' 



C 23 — 1 



S 23 — 



1 TO 3 (™3+TOi)(TO 3 -TO 2 ) 

c 2 3 c (to 3 -toi)(to 3 +to 2 ) ' 



(99) 



r z — 1 - t 2 - 
c 12 — 1 b l2 ~ 



sj 3 to 3 (to|-to^) 
c i3 ^(m?,— to 2 ) 



The combination of the last two equations yields 

sin 2 2#i2 



^ TOito 2 (mj-m{)(mi.-mj) 
(mi+TO 2 ) 2 (to 3 — m2-TO 2 +mim 2 ) 2 



sin 2 2#23 = 4 



(to-i+to 3 -to 2 )(to 2 -toi)(to 3 -to 2 )(to 3 +toi) 
(m 2 -TO 2 -TO 2 +m 1 m 2 ) 2 



(100) 



5.1.1 Numerical analysis 

It is clear from Eqs. (89), (90), and (97)-(100) that this model predicts sf 3 > 0, implying some deviation 
from the TBM limit. The allowed range of values for lepton mixing depends on the value, or range of values, 
used for S13. To perform a numerical analysis let me introduce the parameters k and x defined as 



m 3 - m 2 m 2 
k = , x = — > 1 . 



flOll 



mi 



mi 



One may write all the matrix elements of Vpmns i n terms of these two parameters. In particular s 2 3 may 
be expressed as 
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'13 



x(x — 1) 
(k + l)(k + x- l)(fc + 2x) 



(102) 



The last equation may be used now to invert x in terms of s 2 3 and k, and thus sin 2 2#i2 and sin 2 2#23 may 
be written in terms of s 2 3 and k. A numerical analysis shows that for the range of values 0.033 < s? 3 0.04 
one may obtain large mixing angles for #12 and #23 within the allowed limits of Eq.(l). This region for s 2 3 
is consistent with the upper bound provided by the CHOOZ experiment [T2] (sf 3 < 0.04). I point out the 
allowed magnitudes for lepton mixing in the following (s 2 3 ,fc) parameter space regions: 



global parameter space: 

This region yields the range of mixing angles: 



0.033 s 2 3 ^ 0.04 , 3.2 ^ k sC 4.1 



(103) 



0.64865 ^ sin 2 26 12 ^ 0.818112 , 0.813749 ^ sin 2 29 23 ^ 0.919788 
and the Vpmns unitary mixing matrix with the following range of magnitudes: 



(104) 



VpMNS 



0.830485 
0.254605 
0.401143 



0.874368 
0.371263 
0.432288 



0.442132 
0.765592 
0.367937 



0.526588 
0.773352 
0.461950 



0.181659 
0.524249 
0.784821 



0.2 

0.586563 
0.831963 



(105) 



Recall that the above range of values is restricted by the constraints imposed by the unitarity of Vpmns'-, 
that is, choosing a specific value of one entry further restricts the range of values for the other entries. It is 
clear from Eq.(104) that only part of the values in Eq.(105) are within the allowed limits of Eq.(l). Given 
a particular value for s 2 3 in Eq.(103), it is possible to specify the k parameter region where lepton mixing 
lies within these allowed limits. I point out below these range of values for sf 3 = 0.034 , 0.037 , and 0.04, 
respectively: 



Case A: 



'13 



0.034 , 3.88182 < k < 4.02591 (4.50978 < x < 4.79497) 



0.7 ^ sin 2 2(9 12 < 0.719315 , 0.87 ^ sin 2 2(9 23 ^ 0.878086 



, (106) 



0.859588 - 0.864610 0.467386 -0.476558 0.184390 
Vpmns ~ [ 0.298043 - 0.308732 0.771902 - 0.772539 0.555744 - 0.560673 
0.404500 - 0.407176 0.420784 - 0.429807 0.807246 - 0.810647 



(107) 



Case B: 



'13 



0.037 



3.52059 < k < 3.8732 (4.18727 < x < 4.91525) 



0.7 ^ sin 2 20i2 < 0.749742 , 0.87 < sin 2 2fl 23 < 0.890806 



(108) 



0.849926 - 0.863266 0.466660 -0.490536 0.192353 
Vpmns ~ [ 0.289950 - 0.318087 0.768718 - 0.770414 0.554881 - 0.567795 
0.413159 - 0.420055 0.410422 - 0.434385 0.800381 - 0.809387 



(109) 
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Case C: sj 3 = 0.04 , 3.21323 < k < 3.73671 (3.91261 < x < 5.03901) 

, (110) 

0.7 < sin 2 20i2 < 0.777209 , 0.87 < sin 2 26> 23 0.902305 

0.840573 - 0.861920 0.465933 - 0.503425 0.2 
\ 'rM X s - | 0.282093 - 0.326762 0.765697 - 0.768410 0.554016 - 0.57443 | ( i I I i 

0.421328 - 0.432045 0.400339 - 0.438694 0.793744 - 0.808125 

An additional analysis shows that for 0.0375 < sf 3 < 0.04 and 0.035 < s? 3 < 0.04 one may specify a k 
parameter region where lepton mixing lies within the 3a allowed ranges reported in Refs. [i~3lll4j . respectively. 

5.1.2 Neutrino masses 

With the purpose to obtain some rough estimation for the order of magnitudes of neutrino masses let me 
use the range of values for lepton mixing in Eqs. (106)-(111) and the bounds for Am 2 ol and Am 2 tm of Eq. 
(1). One gets the following neutrino masses. 

Am sol = m-2 — mf = (x~ — 1) my. 

mi » ( 1.796 - 2.145 , 1.750 - 2.320 , 1.706 - 2.494 ) x 10" 3 eV 

m 2 ~ (8.103 - 10.28 , 7.331 - 11.40 , 6.675 - 12.56 ) x 10~ 3 eV , (112) 

m 3 » ( 1.507- 1.892 , 1.349 - 2.039 , 1.215 - 2.188 ) x 10" 2 eV 

where the first, second, and third range of values for each mi ,i = 1,2,3 correspond to sf 3 = 0.034, 
s 2 3 = 0.037, and s 2 3 = 0.04 respectively. 

Am 2 tm = m\ — m\ = k{k + 2x) m\: 

m 1 ss (5.053 - 8.117 , 5.135 - 8.876 , 5.207 - 9.645 ) x 10~ 3 eV 

m 2 ~ (2.279 - 3.892 , 2.150 - 4.363 , 2.037 - 4.860 ) x 10~ 2 eV (113) 

m 3 w (4.240 - 7.160 , 3.958 - 7.801 , 3.710 - 8.464 ) x 10~ 2 eV 



6 Quantitative analysis of quark masses and Vckm 

To leading order in the radiative loop corrections, one gets the approximations 




(114) 
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and then one obtains the relations 



(7 



d 



Y q Y 

1 1 O 1 



_ £12 J i2 ^7>d ^ m d 



16tt 2 



m. 



Y 23 Y 23 rpd m s 

" " i 22 ~ 



1 (irr~ 



<7" = 



to? ~~ 



£l2 2 12 Pli 

16tt 2 1 cr 



r 23 r 23 Z7M 



m t 



(115) 



where the functions , , F^ u , ^r'" are defined analogous to those for the charged lepton sector in 
Eq.(70). 

Hence the mixing angles for the d and u quark sectors, V* and V£, Eq.(64), may be approximated as 



Q d Y 3% F 32 m s 

r>d Y 12 ^23 Q d 

6 12 ~ W W 6 23 

z 23 r 22 



C U r^, £33 £32 ™c 
r 23 ^22 m * 



y9 rpu 

C U, ^ £l2 ^23 n« 

' 12 ~ *23 F 22 23 



(116) 



*33 ™d G d , (m s \2 sf 2 



'13 ~ m s °23 ^ \m b ) s 



yu 



33 m u u 



1 °13 ~ Y u m 

23 *12 m c 



„M I /TOc\2 £12 



6.1 Numerical analysis 

To explore the allowed magnitudes for mixing angles in Vckm and without lost of generality, let me assume 
for simplicity the relationships 



jpd jpu — -77 rpd jpu — -77 jpd jpu — 

^22 — b 22 = J~22 , ^23 ~ ^23 = -^23 , ^32 ~ ^32 = J~32 



F u = T 

(7 — <J 



(117) 



From Eqs. (114)-(117) one gets the following useful relationships to hold 



Y 2 % Y 2 % T 22 



m d m h 
ml 



yl yu t- 
1 12 I 12 £v 

*23 Y 2 U 3 F22 



m u m t 
ml 



yu 
1 12 



m s m u 
m d m c 



yu 
1 23 
yd 
1 23 



m b m c 
m s m t 



(118) 



'12 

b 23 



s 12 r*~> ^12 -^23 

~ Y 2 % T22 



a \2 
b 12 



u yu 
s 23 r-~, 1 33 
-R 

*23 



yd 
1 33 



The combination of Eqs. (114)-(118) yields 



'13 



•yu 

r 33 \ d , 
V d 1 S 13 + 

r 33. 



m t 



^33 A ( m S\2 



'12 



'23 



(119) 



Imposing now for the sake of simplicity 



~yu 
r 33 
•yd 
1 33 



r 23 
yd 
1 23 



= r 



(120) 
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one reaches the simplified relationships between mixing angles in the u and d quark sectors: 




s 23 ~ r s 23 




(121) 



Equation (121) allows one to write the Vqkm = (V£ ) T Vf' quark mixing matrix in terms of four param- 
eters: r and the three mixing angles sf 2 , s^, and sf 3 . A numerical analysis shows that setting for instance 
r = .317239712, corresponding to using the central values for the quark masses m s , nib, m c and mt reported 
in the Particle Data Group, Ref.[K], and the values sf 2 = 0.32721, s^ 3 = 0.0604208 and sf 3 = 0.00921978 
yields the quark mixing matrix (ignoring CP violation): 



Notice that except the matrix element Vtd, the other eight entries lie within the best fit range values reported 
in Ref.[15]. These results suggest that the approach given in Eq.(64) for the orthogonal mixing matrices of 
charged fermions is a good approximation. 

6.2 Quark-Lepton complementarity relations 

Using the quark mixing angles of Eq.(122) and the range of lepton mixing angles of Eqs. (106)-(111) allows 
one to obtain the following rough estimation for the quark- lepton complementary relations [16j : 




(122) 



,PMNS I aCKM 
12 t7 12 



41.543° - 42.152° 



41.543° - 43.139 



41.543° - 44.066° 



(123) 



\PMNS I oCKM 
23 a 23 



36.835° - 37.184' 



O 



36.835° - 37.754' 



O 



36.835° - 38.295 



o 



for s? 



13 — 



0.034 , 0.037 , and 0.04, respectively. 
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7 FCNCs and rare decays for charged leptons 



The new exotic scalar particles introduced to implement the radiative mass generation mechanism have the 
capability to induce FCNCs and contribute to "flavor violation" processes such as F — > /1/2/3, to "radiative 
flavor violating" processes such as \i — * ej, r — * fij, r — > ej, as well as to the "anomalous magnetic 
moments" (AMMs) of fermions. In this section I compute roughly these additional contributions for the 
charged leptons. 

Once the generation of fermion masses is completed, the transformations between gauge (0 superscript) 
and mass (physical) eigenstates are for scalars <i>j = UijCFj, Eq.(17), for charged leptons 

^°eL,eR = V e L,eR 1peL,eR , (124) 

where 

VeL,eR = V^V^eR > ^eR = (e° ', T°) L ,R , ^ = M, t) L ,R , (125) 

and analogous transformations for quarks. 

7.1 Lepton flavor violation (LFV) processes F — ► /1/2/3 

The scalar fields (<pg, 0io, 4>i2, <t>il) allow tree level flavor changing vertices through the couplings in Eq.(7). In 
particular they may induce tree level "lepton flavor violation" (LFV) processes such as r — > fifi/j,, r — ► /i/ie, 
t — > fj.ee, t — > eee, and [i — * eee. The generic diagram for these processes is shown in Fig. 4. The decay 
rate contribution from this generic diagram may be taken as [17] 

m 5 Y 4 

4> 

with Y\ being a coupling constant. 



7.1.1 (j, — * eee 

Here I discuss some details about the decay \i — ► eee. This rare decay is of particular interest to be analyzed 
because experimentally it is strongly suppressed. The dominant contribution to this decay comes from the 
diagrams of Fig. 5. Then, from Eqs.(7) and (126), a rough estimation for this decay rate may be written as 

and therefore the branching ratio for this process is^l 

BR^ -> eee) = - l£{ }, eee < 1 x 1(T 12 , (128) 

where I take 

- K^) 4 1^? <i29) 

3 I write the experimental bounds reported in Particle Data Group Ref.[15]. 
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7.1.2 r — * fififi 



(130) 



Using the mean life of r = (290.6 ± 1.0) X 10~ 15 s, and hence (IY) T = ± ps 2.2711631 x 10~ 12 GeF, one gets 
the branching ratio 



5i?(r MMM) ~ am r 4 y 2 4 3 { } T/W < 1.9 x 10-' , 



(131) 



with C T = 4.107105 x 10 6 



7.1.3 r — > [i + [Te 



T(r — > e 



m5(yi 2 y 23 ) 2 1 

3072 vr 3 2 



(132) 



with the branching ratio 



BR(t~ m + M~0 ps C T m T 4 (y 12 y23) 2 { } W e < 2 x lO" 7 



(133) 



7.1.4 r — > ^ + e e 



T(r — > /i + e e 



^(yi 2 y 23 ) 2 1 

3072 vr 3 2 



(^(e^) 2 +( ^(e^) : 



(134) 



rfiee 



BR(j- /x+e-e") « am r 4 (yi 2 y 23 ) 2 { } TMee < 1.1 x 10~ 7 



(135) 



7.1.5 r — ► eee 



r(r — ► eee) 



3072 vr 3 2 | [\< ~Mf ) [1. 



(136) 



BR(t eee) ps C T m T 4 (yi 2 y 23 ) 2 { } reee < 2 x 10" 



(137) 



7.2 Anomalous magnetic moments and radiative rare decays F — > f-y 

The amplitude for the radiative process f\ — ► / 2 7 with /i and / 2 being two equally charged fermions and 7 
a real photon is written as [18] 



iM(h( Pl ) - / 2 (p 2 ) + 7) = »«2(P2) eVl ( ) 5 /i/2 + 



^ **7nU. . ^ gAgj^. (F V (0)+F A (Q) )) //|(/)|) (|3 , S) 
m\ + m 2 
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where F 2 gives the AMM (electric dipole moment) for the fermion f\ when f\ = f 2 - The generic diagrams 
for the process fn — > f 2 R 7 are shown in Fig. 6, in these diagrams a stands for a mass eigenstate scalar 
field. The respective evaluation of these diagrams gives 

y2 Y 2 a 
iA L & ■^N(M k ,m i )e 2 Ri<Ti a 'q u e ll eiL and iA R » -^N(M k , mi)e 2L i^ v q^em , (139) 

where the second amplitude comes from the diagrams where L and R are interchanged, and N(M k ,rrii) may 
be approximated as 



rrii_. Mi 



N{M k , mi )^^\n^- (140) 
Mr mf 



in the limit M k 3> m*. Notice that due to scalar field mixing the contribution of these loops is finite as 
those in the mass case. 

Because of the fermion mixing matrices structure the diagrams that make the largest contribution to 
the AMMs of the charged leptons are, for the electron, the diagram with the muon inside the loop, and for 
the muon and tau, the diagrams with tau as the internal fermion. 

7.2.1 Muon anomalous magnetic moment 

The dominant contribution for the muon AMM comes from the diagram of Fig. 7, where the insertion of 
a photon on the internal lines is understood as in the generic diagrams of Fig. 6. The expression for this 
scalar contribution is [I8l [19] 



a, = ^(V eL ) 22 (V eR ) 22 (G^ + G^)^^(G^ + G^ 

^ m^m T Yj 3 ^ U lk U 4k , Mj 



(141) 



where 



qui = G m = U lk U 4k (V eL ) 3l (V eR ) 3l N(M k , mi ) « Z k U lk U 4k N(M k ,m 7 

~ m r l^ k ~^r~ in -^r . 

7.2.2 Electron and tau anomalous magnetic moments 

Performing a similar analysis for e and r leptons, one gets 
For electron: 

f^eXvi ,„ e , , ^c Bl . rn e m^ Y? 2 ^ U 2k U 3k , M| 



(142) 



16vr 2 v ; 8tt 2 Ml m\ 



where 
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G e L = Q e R = ^ U 2k U 3k (V eL ) 2i (V eR ) 2i N(M k , mi ) ~ J2 k U2kU 3k N(M k ,m ll ) 



(144) 



For tau: 



(G TL + G-) « ^ V ™* In ^ (145) 
16vr 2 v 7 8vr 2 ^ M 2 ml y 1 

k K 



where 



G r L = G r R = U 2k U 3k (V eL ) 3i (V eR ) 3i N(M k , mi ) « £ U 2k U 3k N(M k , m T ) 

~ m T > — 779— In — f- . (146) 

7.2.3 Radiative decay /x — > e7 

A similar analysis to the one for the muon AMM leads to the decay rate 

r(M _ e7 ) = K + ^ y/ fl - ^) 2 (1 - ^) (|(V eL ) 22 (^)nG™| 2 + |(Vk)n(V^) 22 G^| 2 ) 
- igffl 2 (|G^| 2 + |G^| 2 ) 



(16) :1 7T b 



G™ = Z k , i Uv e U M (V eL )si(VeR)2iN(M k ,mi) ~ (^23 m T E fc ^^Mn^l 
G"*«* = Y, kt U 2k U 4k (V eL ) 2t (V eR ) 3i N(M k ,m i ) « (V eL ) 23 m T J2 k ^^\n^ 

k r 

The resulting branching ratio may be expressed as 

3 m 2 /il%\ 4 



(147) 



(148) 



BR(ji - e 7 ) « — — ^ (Y 12 Y 23 Y{ } Me7 < 1.2 x lO" 11 . (149) 



2vr 2 mf l \g w J 
7.2.4 Radiative decays r — > ^7 and r — > e 7 

Carrying out a similar analysis, one gets 

T -► ^7: 
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r (^^) « m (igy (iG^^i 2 +|g^i 2 ) 

f2 \ 2 / rr r , ,,2 \ 2 



w 

(16)3 ttS 1 ^fc M 2 m + ^fc Mf 111 



where 



2 

M, 2 



and branching ratio 

5i?(r ^7) ~ C>r (^23^33) 2 { }r M7 < 6.8 X 1(T 8 , 

with CL = 6.24205152 x 10 5 . 



t — > e^f 

(16) a w 



r(r-e 7 ) « m (S^ 3) (|G^ e *| 2 + |G^| 2 ) 



where 



G TL£R = Ek i Ui k U 3k (V eL ) 2l (V eR ) 2t N(M k , rm) « m M £ fc In $ 

k L 1 

G TReL = Ek,i U 2k U 4k (V eL ) 2i (V eR ) 2i N(M k , rm) « m M £ fc ^ ln 5 



and branching ratio 



BR(t -► e 7 ) « C> 2 m 2 (l^s) 2 ! }re 7 < 1.1 x 1(T 7 
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8 Summary and conclusions 



I have reported a detailed analysis on fermion masses and mixing, including neutrino mixing, within the 
context of an extension of the standard model with an U(1)h flavor symmetry and hierarchical radiative 
mass mechanism [2] . The results of this analysis show that this model has the capability to accommodate 
the observed spectrum of quark masses and mixing angles in the Vckm, as it is shown through the analysis 
in Sees. 3 and 6. In a similar way the spectrum of charged lepton masses is consistently generated through 
the analysis presented in Sees. 3 and 4. Upper bounds for the charged lepton mixing angles are given 
in Eq. (83). These upper bounds imply that mixing in the lepton sector comes almost completely from 
neutrino mixing; that is, VpMNS ~V U . In this approach all lepton mixing elements in Vpmns are written 
completely in terms of neutrino masses. A numerical analysis shows that using 0.033 < sf 3 < 0.04 one gets 
large mixing angles for 8 12 and <9 23 , 0.7 < sin 2 28 12 < 0.777209 and 0.87 sC sin 2 29 23 < 0.902305, within the 
present allowed 3a limits as reported by recent global analysis of neutrino data oscillation (8J fl~3l fT4"] . Using 
these allowed ranges of values for quark and lepton mixing, predictions for neutrino masses and quarkdepton 
complementary relations are given in the Eqs. (112), (113), and (123), respectively. 

From the phenomenological point of view it is interesting to look for a set of scalar mass parameters in 
M 2 , Eq. (165), that allows us to account for the strong experimental suppression on LFV processes, such 
as [i — > eee, radiative rare decays fi — > ej, r — * ^7, r — ► ej, and the muon anomalous magnetic moment. 
To achieve this goal a detailed numerical analysis and fit it is needed, trying to keep at least the lowest 
scalar mass eigenvalue rj\ within few TeV 2 . However, it is important to comment that Eq. (83) gives a 
good approximation for the upper bounds on charged lepton mixing angles, and hence Vpmns ~ V v would 
remain as a good approach in this model. 

Thus, the contribution of my analysis in comparison to the one realized in Ref. [2] may be summarized 
in the following aspects: 

• Scalar sector: I have performed the analysis by considering the most general structure for M?. 

• Charged fermion sector: 

— I have obtained and then diagonalized the quark and lepton mass matrices at one and two loops 
in close analytical form. 

— Taking advantage of the strong hierarchy of quark and charged lepton masses, approximate 
expressions for the orthogonal mixing matrices of charged fermions are obtained. 

— I have reported general analytical expressions for the branching ratios of LFV processes, radiative 
rare decays and for the AMMs of charged leptons. 

• Neutrinos: The Vpmns lepton mixing matrix is obtained and written completely in terms of the 
neutrino masses, and numerical results for lepton mixing angles are provided. 
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9 Appendix: Diagonalization of a generic real symmetric 3x3 mass 
matrix 



In this appendix I give the details to diagonalize a generic real symmetric mass matrix defined as 



M 




(156) 



One can diagonalize this matrix M through the orthogonal matrix V as V T MV = diag{\\, A2, A3), Aj, 
i = 1,2,3 being the eigenvalues of M. The determinant equation det\M — A| = imposes the constraint that 
each one of the eigenvalues A, obeys the cubic equation 

- A 3 + (a + b + c)A 2 - (ab - d 2 +ac-e 2 +bc- f 2 )X + abc - fa - e 2 b - d 2 c + 2def = . (157) 

Thus, from Eq. (157) one obtains the following nonlinear relationships to hold: 

Ai + A 2 + A 3 = a + b + c 
AiA 2 + A1A3 + A 2 A 3 = ab - d 2 + ac - e 2 + be - f 2 
A1A2A3 = abc- f 2 a - e 2 b - d 2 c + 2def . 



(158) 



I do not impose any hierarchy between the eigenvalues Ai, A2, and A3. However, I assume they are nonde- 
generated. Computing now the eigenvector^], the orthogonal matrix V may be writing as 



V 



X 



J\(A 2 ) F 2 (A 3 ) \ 
y A 2 (A 2 ) z A 3 (A 3 ) 



^i(Ai) 
Ai(Ai) 



y 



z 



F 3 (X 3 ) 
A 3 (A 3 ) 



T W F 3(A 2 ) 
V X Ai(Ai) y A 2 (A 2 ) 



(159) 



where x, y, and z are normalization constants, and the functions involved are defined as 



Ai(A) 


= (6- 


A)(c- 


A) 


-f 


Fi(A) = 


-d(c- 


A)+e/, 


A 2 (A) 


= (0- 


A)(c- 


A) 


-e 2 


F 2 (A) = 


-e(b- 


A) + df , 


A 3 (A) 


= (a- 


A)(6- 


A) 


-d 2 


F 3 (A) = 


-f(a~ 


A) + de . 



(160) 



Using properly Eqs. (157) and (158), it is possibly to check the orthogonality between columns(eigenvectors) 
of V. Moreover, the functions Aj(A) and Fi(X) in Eq.(160) satisfy the important and useful relationships 



4 Here still unnormalized 
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if(A) = A 1 (A)A 2 (A) , F 1 (A)F 2 (A) = A 1 (A)F 3 (A) , 

Ff(A) = Ai(A)A 3 (A) , F 1 (X)F 3 (X) = A 2 (X)F 2 (X) , (161) 
F|(A) = A 2 (A)A 3 (A) , F 2 (A)F 3 (A) = A 3 (A)iq(A) . 



Defining 



h{\) = Ai(A) + A 2 (A) + A 3 (A) = 3A 2 - 2(o + b + c)A + ab-d 2 +ac-e 2 +bc- f 

= 3A 2 - 2(Ai + A 2 + A 3 )A + A X A 2 + AiA 3 + A 2 A 3 , 



(162) 



explicitly, 



MAi) = A 2 -(A 2 + A 3 )A 1 + A 2 A 3 = (A 2 -A 1 )(A 3 -A 1 ), 

h{\ 2 ) = A|-(Ai + A 3 )A 2 + AiA 3 = (Ai-A 2 )(A 3 -A 2 ), (163) 
h(X 3 ) = A|-(Ai + A 2 )A 3 + A 1 A 2 = (A 1 -A 3 )(A 2 -A 3 ). 

One can use now the relationships of Eqs. (161) and (163) to normalize the eigenvectors, obtaining that in 
general the square matrix elements V 2 - i,j = 1, 2, 3 may be expressed as 



^=w-° andhence (164) 

Equation (164) defines the magnitudes for the matrix elements Vij in Eq.(159). Setting now the diagonal 
elements of V as positives, x > 0,y > 0, and z > 0, the signs of the off diagonal elements, Vij, i ^ j, may 
be obtained directly from the Eq.(159) in a particular set of giving parameters a, 6, c, d, e, and / that define 
the real symmetric mass matrix M in Eq. (156). 

It is important to mention here that the method introduced in this Appendix to diagonalize a generic 
3x3 real symmetric mass matrix agrees with the diagonalization performed in Ref. [2D] for the special case 
of Fritzsch's ansatz, a = e = 0. 



9.1 Diagonalization of the generic exotic scalar mass matrices 



The most general square scalar mass matrix for the exotic scalar fields, which mediate the radiative mass 
generation of the light fermions at one and two loops in the u, d, and e charged fermion sectors, may be 
written as 



/ 


a i 


b 





\ 




b 


«2 


c 










c 




d 


V 








d 


o 4 / 



(165) 
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This matrix may be diagonalized through the orthogonal matrix U as U T M^U = cfoci<?fai, 774), = 

Mf ,i = 1, 2, 3, 4 being the eigenvalues of M|. Using the same procedure and method introduced previously, 
the orthogonal matrix U may be writing as 



U 



x i i him) z i him) hM \ 

■' A 2 fa 2 ;• " "- 1 • 



X 



i him) 

Aifai) 



y 



A 3 (r7 3 ) 

/ g3iv3 



Aiirji) 



. , + i .94(774) 

Z A 3 (7? 3 ) 1 A 4 (7 ?4 ) 



X 



1 him) „,/ 9zim) 



Aifai) 



y 



A 2 fa 2 ) 



V X Ai(»ji) y A 2 fa 2 ) Z Agfa) 



,/ h 4 jm) 

J A 4 (?74) 



where x', y', z', and t' are normalization constants, and the functions involved are defined as 



and 





Alfa) 


= (a 2 


-7/) (03 - 


r/) (04 - 


7?) - (a 2 - 


- 7])d 2 - 


(a 4 - 77) c 2 , 




A 2 fa) 


= (ai 


- »/) [(as 


- »/)(04 


-r?)-^] 








A 3 fa) 


= (a4 


-»/) [( a i 


- V)(02 


-r?)-6 2 ] 








A 4 fa) 


= (ai 


-??)(a 2 - 




77) - (ai - 


- 7?)c 2 — 


(a 3 - 7?)6 2 , 


hin) 


= ~b [(03 


- »/)(a4 


-77)-^ 


] > 


53 fa) = 


-c(ai - 


7?)(a 4 -77) , 


hin) 


= 60(04 — 


77) 






54 fa) = 


cd(ai — 


??) . 


hiv) 


= —feed 






5 


Mv) = 


-d [(ai 


- r))(a 2 -77) 



(166) 



(167) 



(168) 



These functions satisfy relationships analogous to those of Eq.(161), allowing us to obtain the normalization 
constants and then to write the square matrix elements as 



2 = A^ 

where 

Km) = im - m)(m - m)im - m) 
Hm) = ivi - m){m - m)im - m) 

and leads to the useful equalities 



and hence 



\Uij\ = 



'Aifaj 



(169) 



Km) = im - m)im - m)(m - m) , 
Km) = im - m)im - m)im - m) , 



(170) 



hiVk 

Kvk 



gffafc) 
Km) 



Mm 
Km 



ffjjfa* 

Kvk 



(171) 
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Figure 1: Generic diagrams contributing to fermion masses, (a) Dirac-type couplings, (b) Majorana-type 
couplings. 
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Figure 2: Mass diagrams for the charged lepton sector. 
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(a) (b) 
Figure 5: Diagrams for the rare decay \i — ► eee. 




Figure 6: Generic diagrams for the process l\ — > Z27, where a scalar mass eigenstate a is involved. 
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M-2 x2 ^ " K e iL e iR K 

Figure 7: Main contribution to the muon anomalous magnetic moment (the insertion of a photon on the 
internal lines is understood as in the Fig. 6). 
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